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Abstract. In this article we show how to compute the semi-classical spectral 
measure associated with the Schrodinger operator on R", and, by examining 
the first few terms in the asymptotic expansion of this measure, obtain inverse 
spectral results in one and two dimensions. (In particular wc show that for 
the Schrodinger operator on with a radially symmetric electric potential, 
V, and magnetic potential, B, both V and B are spectrally determined.) We 
also show that in one dimension there is a very simple explicit identity relating 
the spectral measure of the Schrodinger operator with its Birkhoff canonical 
form. 



be the semi-classical Schrodinger operator with potential function, V{x) G C°°(R"), 
where A is the Laplacian operator on R". We will assume that V is nonnegative 
and that for some a > 0, V^^{[0,a]) is compact. By Friedrich's theorem these 
assumptions imply that the spectrum of Sn on the interval [0, o) consists of a finite 
number of discrete eigenvalues 



with N{h) oo as h ^ 0. We will show that for / e C°°(M), with supp(/) C 
(— oo,a), one has an asymptotic expansion 



1. Introduction 



Let 



(1.1) 



Sn = --A + V{x) 



(1.2) 



Ai(/i), 1 < i < N{h) 



OO 



(1.3) 




with principal term 



(1.4) 
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and subprincipal term 



(1-5) ^i(^) = "24 / ^^'^(y + ^^""^^ E ^ 

i * 

We will also give an algorithm for computing the higher order terms and will show 
that the /c*'' term is given by an expression of the form 

(1.6) Mf)= [ E f^^'HY + ^i'^^^PkADV,--- ,D''^V)dxdC 

where pkj are universal polynomials, and D'^V the A;*'' partial derivatives of V. 
(We will illustrate in an appendix how this algorithm works by computing a few of 
these terms in the one-dimensional case.) 

One way to think about the result above is to view the left hand side of (1.3) 
as defining a measure, fin, on the interval [0, a), and the right hand side as an 
asymptotic expansion of this spectral measure as h ~f 0, 

(1-7) ^'^-T.f'""(iT 



where Hk is a measure on [0, a) whose singular support is the set of critical value of 
the function, V. This "semi-classical" spectral theorem is a special case of a semi- 
classical spectral theorem for elliptic operators which we will describe in §2, and in 
§3 we will derive the formulas (1.4 1 and (1.51 and the algorithm for computing (1.6) 
from this more general result. More explicitly, we'll show that this more general 
result gives, more or less immediately, an expansion similar to (1.7 1, but with a 
"(d_-^ikv pig^pg Qf ^Yie {j^f''■ We'll then show how to deduce (O) from this 
expansion by judicious integrations by parts. 

In one dimension our results are closely related to recent results of |Col05j . 
|Col08| . |CoG] . and |Hezj . In particular, the main result of |CoG| asserts that if 
c € [0,a) is an isolated critical value of V and V~^{a) is a single non-degenerate 
critical point, p, then the first two terms in (1.7) determine the Taylor series of V 
at p, and hence, if V is analytic in a neighborhood of p, determine V itself in this 
neighborhood of p. In ^ColOS] Colin de Verdiere proves a number of much stronger 
variants of this result (modulo stronger hypotheses on V). In particular, he shows 
that for a single-well potential the spectrum of Sn determines V up to V{x) <-> 
V{—x) without any analyticity assumptions provided one makes certain asymmetry 
assumption on V . His proof is based on a close examination of the principal and 
subprincipal terms in the "Bohr-Sommerfeld rules to all orders" formula that he 
derives in |Col05j. However, we'll show in §4 that this result is also easily deducible 
from the one-dimensional versions of (1.4) and (1.5), and as a second application 
of (1.4) and (1.5), we will prove in §5 an inverse result for symmetric double well 
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potentials. We will also show (by slightly generalizing a counter-example of Colin) 
that if one drops his asymmetry assumptions one can construct uncountable sets, 
{Va,ct G (0, 1)}, of single- well potentials, the Vq's all distinct, for which the /i^'s in 
(1.7l are the same, i.e. which are isospectral modulo 0{1t"). 



In one dimension one can also interpret the expansion (1.7) from a somewhat 
different perspective. In §6 we will prove the following "quantum Birkhoff canonical 
form" theorem: 

Theorem 6.1. If V is a simple single-well potential on the interval V^^ ([0,a)) 
then on this interval Sfi is unitarily equivalent to an operator of the form 

(1.8) HQB{St\h') + Oih°^) 

where S'^"'^ is the semi-classical harmonic oscillator: the 1-D Schrodinger operator 
with potential, V{x) — ^ . 

Then in §7 we will show that the spectral measure, fi^, on the interval, (0, a), is 
given by 

(1.9) Mf) = £ fit)^it,f^')dt 
where 

(1.10) HQB{s,h^) = t ^ s = K{t,h^). 

In other words, the spectral measure determines the Birkhoff canonical forms and 
vice-versa. 



The last part of this paper is devoted to studying analogues of results (1.3l-(1.7) 
in the presence of magnetic field. In this case the Schrodinger operator becomes 

(1.11) 4'")=^(^A+„,(,))2 + y(,) 

^-^ I axk 



fe=i 



where a = ^ audxk is the vector potential associated with the magnetic field and 
the field itself is the two form 

(1.12) B ^ da^ ^ Bi-jdxi A dx^. 

For the operator the analogues of ( |1.3| -( 1.7) are still true, although the 

formula (1.6) becomes considerately more complicated. We will show that the 



subprincipal term (1.5) is now given by 

(1-13) ^ / f^^\\Y.^i^ + a,f + V{x)){^2Y,^^^+\\Br)dxdt 

As a result, we will show in dimension 2 that if V and B are radially symmetric 
they are spectrally determined. 
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2. Semi-classical Trace Formula 

Let 

(2.1) Pn= Yl a^{^M^DX 

|Q|<r 

be a semi-classical differential operator on M" , where (x, fi.) G C°° (M" x M) . Recall 
that the Kohn-Nirenberg symbol of Pn is 

(2.2) p(a;,e,/i) = ^a„(x,;i)r 

a 

and its Weyl symbol is 

ih 

(2.3) p'"(x,e,/i) =cxp(--D^9,)p(a:,^,S). 

We assume that is a real-valued function, so that Pri is self-adjoint. Moreover 
we assume that for the interval [a, 6], (p"')^^([a, 6]), Q <h< ho, is compact. Then 
by Friedrich's theorem, the spectrum of Pn, h < Hq, on the interval [a,b], consists 
of a finite number of eigenvalues, 

(2.4) \,{H), l<i< N{h), 
with N{h) ^ oo as /i — > 0. Let 

(2.5) p{x,O^p{x,^,0)^p^{x,^,0), 

be the principal symbols of Pfi. 

Suppose / £ C^(K) is smooth and compactly supported on (a, 6). Then 

f{Pn) = -L / /(t)e"^« dt, 
V 27r J 

where / is the Fourier transform of /. 

Theorem 2.1 ( [GuS] ). The operator f{Pn) is a semi- classical Fourier integral 
operator. In the case p{x, ^, h) =p(x,S,), i-C. aa{x,h) are independent of h, f{Ph) 
has the left Kohn-Nirenberg symbol 

bf{x, i,h)^j2 ( E ^) ((7 ^)'^) ^)) I ■ 

k \l<2k J 



(2.6) 
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It follows that 
(2.7) 



trace/(Pft) = h''' J bf{x,^,h) dxd^ + 0{h°°). 



The coefficients bk^iix,^) in (2.6 1 can be computed as follows: Let Qa be the 
operator 



(2.8) 



\ [dx + it ^ 
a\ \ ox 



Let bk{x,^,t) be defined iteratively by means of the equation 
ldb„ 



(2.9) 

with initial conditions 

(2.10) 

and 

(2.11) 



i dt 



\a\>l k+\a\=m 



6„(x,C,0) =0 



for TO > 1. Then it is easy to see that bk{x,£,,t) is a polynomial in t of degree 2k. 
The functions bk^i{x,S,) are just the coefficients of this polynomial, 



(2.12) 



bk{x,^,t)^ J2bkAx:0t'- 



l<2k 



3. Spectral Invariants for Schrodinger Operators 
Now let's compute trace/(S'fi) for / S C^{—a,a) via the semi-classical trace 



formula (2.7). Notice that from (2.6 1, (2.7 1 and (2.10) it follows that the ffist trace 



invariant is 



/(p(a:,0) dxd^, 

which implies Weyl's law, f |GuS] §9.8), for the asymptotic distributions of the 



eigenvalues (2.4) 



To compute the next trace invariant, we notice that for the Schrodinger operator 



(l.ll, 



(3.1) 

so the operator Qa becomes 
(3.2) 



1 dW 
ox 
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It follows from that 

|q!|>1 fc+|a|=m 



E 



i \dxk dxk 



brn. — 



k ^ ^ k 

Since 6o(x,f,t) = 1 and &i(x,^,0) = 0, we have 



— it—\ 
dxk dxk ) 



dV 

dxi ' 



and thus 
i dt 



E 



^k f d 



dV ^ 1 



yxk d^kj 2 ^ UXi 



dxk dxk ) 



(1) 



y E ^'^^^ 



k.l 



it- 



dV dV 



dxkdxi dxk dxi 



d^V ,^ dV dV 

dxk dxk 



dxi 



It follows that 
(3.3) 

^^(-'^'*)-4Ea4+T 

k 



8 ^ ' dxk dxi ' 



Thus the next trace invariant will be the integral 
(3.4) 



^)) 



6 ^ ^'^^ dxhdx 

k,l *- 



We can apply to these expressions the integration by parts formula, 



(3.5) 
and 



A{x)^B'i^ + V{x)) dxd^ 
dxk 2 



(3.6) j ^kiiA{x)B'{^^ + V{x)) dxd^ = ~J 5iA{x)B{^^ + V{x)) dxd^. 



Applying (3.5 1 to the first term in (3.4 1 we get 

C2 
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and by applying (3.6 1 the fourth term in (3.4) becomes 



J s^^dxj ^ ^ 



e 



V{x)) dxdC 



Finally applying both (3.6 1 and (3.5 1 the third term in (3.4| becomes 

C2 



So the integral (3.4 1 can be simplified to 



2i/i:(£)=/"(| + i'W)<i.<ie 



We conclude 



Theorem 3.1. The first two terms of (2.1) are 
(3.7) 



traeef{Sn) = J f{j+V{x))dxd^+^h'J Y,{^f fi^){^l+V{x)) dxdi+0{h^). 



In deriving (3.7 1 we have assumed that / is compactly supported. However, 
since the spectrum of Sn is bounded from below by zero the left and right hand 
sides of ( |3.7[ ) are unchanged if we replace the "/" in (3.7 1 by any function, /, with 
support on (— oo, a), and, as a consequence of this remark, it is easy to see that the 
following two integrals. 



(3.8) 

and 
(3.9) 



+ V{x)<X 



dxd^ 



L 



y^i^fdxd^ 



are spectrally determined by the spectrum (2.4 1 on the interval [0, a]. Moreover, 
from (3.7), one reads off the Weyl law: For < A < a, 

(3.10) #{A,(fi) < A} = (27rfi)-" (vol(^ + V{x) < A) + 0{h) 



We also note that the second term in the formula (3.7) can, by (3.6), be written 
in the form 



24 

and from this one can deduce an fi^-order "cumulative shift to the left" correction 
to the Weyl law. 
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3.1. Proof of (1.6). To prove (1.6 1, we notice that for m even, the lowest degree 
term in the polynomial bm is of degree ^ + 1, thus we can write 



m+l 



Putting this into the the iteration formula, we will get 



m + l, -^^kdbm-i.i 
— -. — o,n,i = 



dxk 



-y 

2 ^ 



Dm-2,i+l 

dxl 



_ i d dV 
2 dxk dxk 



dV d 



)bm- 



dxk dxk 

from which one can easily conclude that for / > 0, 



2,1 ■ 



5E( 



dV . 
dxk' 



2,1-1, 



(3.11) 



'Jm,l 



Er(|;)v.(w,. 



,D"'V) 



where Po,,/? is a polynomial, and \a\ 



> 2/ — 1. It follows that, by applying the 



integration by parts formula ( 3.5 1 and ( 3.6 1 , all the /(™+') , ^ > 0, in the integrand of 
the fi."th term in the expansion (2.6 1 can be replaced by /(™). In other words, only 



derivatives of / of degree < 2k figure in the expression for Vk{f)- For those terms 
involving derivatives of order less than 2k, one can also use integration by parts to 
show that each J*^™^ can be replaced by a and a In particular, we 



can replace all the odd derivatives by even derivatives. This proves (1.6 1 



4. Inverse Spectral Result: Recovering the Potential Well 

Suppose y is a "potential well", i.e. has a unique nondegenerate critical point 
at a; = with minimal value V^(0) = 0, and that V is increasing for x positive, and 
decreasing for x negative. For simplicity assume in addition that 

(4.1) - V'{-x) > V'{x) 



and (3.91 to 



holds for all x. We will show how to use the spectral invariants (3 
recover the potential function V{x) on the interval |a;| < a. 

For < A < a we let ~X2{\) < < a;i(A) be the intersection of the curve 

%2 + y{x) = A with the x-axis on the x — ^ plane. We will denote by Ai the region 
in the first quadrant bounded by this curve, and by A2 the region in the second 



quadrant bounded by this curve. Then from (3.8) and (3.9) we can determine 



(4.2) 

and 
(4.3) 



Ax J A2 



dxd^ 



f + f V'{xfdxd£,. 

J Ai J A2 
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V = V{x) 




+ V{x) ^ A 



Figure 1 . Single Well Potential 



Let X = /i(s) be the inverse function oi s — V{x)tX G (0,a). Then 



V'{xf dxdS, = / V'{xY / didx 

Ai Jo Jo 



xi{\} 



V'{xfy^2X-2Vix) dx 



Similarly 



V2X-2sV'{fi{s)) ds 



V'{xf dxdS, = / V2A - 2s ( ^ ) ds, 



IA2 JO 

where x — /2(s) is the inverse function of s = V{—x),x G (0,a). So the spectrum 
of Sfi determines 

(4.4) 



VA^((?^r + (?)-M ds. 



ds 



ds 



Similarly the knowledge of the integral (4.2 1 amounts to the knowledge of 



(4.5) 



+ ds. 
ds ds 



Recall now that the fractional integration operation of Abel, 



(4.6) 



1 



r(a) 



iX-tr~'git) dt 



for a > satisfies J^J** = J''^^. Hence if we apply J^/^ to the expression (4.5 1 and 
(4.4 1 and then differentiate by A two times we recover ^ + ^ and (^)^^ + (^)~^ 
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from the spectral data. In other words, we can determine /[ and /2 up to the 
ambiguity /{ ^ /2- 

However, by (4.1 1, /( > /j. So we can from the above determine /{ and /j, and 
hence fi,i — 1,2. So we conclude 

Theorem 4.1. Suppose the potential function V is a potential well, then the semi- 
classical spectrum of Sn modulo o{h^) determines V near up to V{x) <-> V{—x). 



Remark 4.2. The hypothesis (4.11 or some "asymmetry" condition similar to it is 
necessary for the theorem above to be true. To see this we note that since V{x) and 
V{—x) have the same spectrum the integrals in (1.6 1 have to be invariant under 
the involution, x 
Now let y : M ^ 
X —^ ±cxD and 



(a) 



and 



(6) 



~x. (This is also easy to see directly from the algorithm ( |2.9[ ).) 
be a single well potential satisfying V{0) ~ 0, V{x) — > +oo as 

V{-x) = V{x) for fc > and 2k < x < 2k + 1 
V{-x) < V{x) for A: > and 2fc + 1 < x < 2A: + 2. 



Now write the integral ( |1.6[ ) as a sum 
(4-7) E 



where Ik is the set, {{x,^), k < x < k + 1}, and for a E (0, 1) having the binary 
expansion .aia2a3 • • • , = or 1, let Vq, be the potential 

Va{x) = V{x) on 2k < X < 2k + Hi ak ^ 



and 



Vaix) = V{-x) on 2fc < x < 2fc + 1 if flfc = 1. 



In view of the remark above the summations (4.7 1 are unchanged if we replace V 
by Vc.. 



Remark 4.3. The formula (4.5 1 can be used to construct lots of ZoU potentials, i.e. 
potentials for which the Hamiltonian flow vh associated with H ~ + V{x) is 
periodic of period 2tt. It's clear that the potential V{x) = x"^ has this property and 
is the only even potential with this property. However, by (4.5) and the area-period 
relation (See Proposition 6.1) every single-well potential V for which 

fi{s) + f2{s)^2s^'^ 

has this property. We will discuss some implications of this in a sequel to this 
paper. 
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5. Inverse Spectral Result: Recovering Symmetric Double Well 

Potential 

We can also use the spectral invariants above to recover double-well potentials. 
Explicitly, suppose y is a symmetric double-well potential, V{x) = V{—x), as 
shown in the below graph. Then V is defined by two functions Vi, V2: 




Figure 2. Double Well Potential 



As in the single well potential case, let /i, /2 be the inverse functions 
and 

X = f2{s) s = 1/2(2; - a). 

For A small, the region {(a;,C) I 2 ^(■^) — ^} ^ indicated in figure 2, so if 
we apply the same argument as in the previous section, we recover from the area 
of this region the sum ^ + ^ via Abel's integral. Similarly from the spectral 
invariant Ja^ ^yi^^-^^x^^'^"^ dxd^ we recover the sum (^)~^ + (^)~^- As a result. 



we can determine Vi and V2 modulo an asymmetry condition such as (4.1). 

The same idea also shows that if V is decreasing on (—00, —a) and is increasing 
on (6,00), and that V is known on (~a,b), then we can recover V everywhere. In 
particular, we can weaken the symmetry condition on double well potentials: if V 
is a double well potential, and is symmetric on the interval y^^[0, V^(0)], then we 
can recover V. 



6. The Birkhoff Canonical Form Theorem for the 1-D Schrodinger 

Operator 

Suppose that ([0, a]) is a closed interval, [c, d], with c < < d and V{0) = 0. 
Moreover suppose that on this interval, V" > 0. We will show below that there 
exists a semi-classical Fourier integral operator, 

U : C^'{M.) C°°(M) 
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with the properties 

(6.1) Uf{Sn)U' = fiHQsi.Sl'^^h)) + 0{fi^) 
for all f €Cg° ((-oo, a)) and 

(6.2) UU*A = A 

for all semi-classical pseudodifferential operators with microsupport on ([0, a)). 
To prove these assertions we will need some standard facts about Hamiltonian 

systems in two dimensions: With H{x, C) = V + ^(^) ^ above, let v = Vh he the 
Hamiltonian vector field 

_dH d dH d 
dx dx 

and for A < a let 7(4, A) be the integral curve of v with initial point, 7(0, A), lying 
on the X-axis and if(7(0. A)) = A. Then, since L^H = 0, H{'y{t,X)) = A for all t. 
Let T{X) be the time required for this curve to return to its initial point, i.e. 

7(t, A) 7^ 7(0, A), for < t < T(A) 

and 

7(T(A),A) = 7(0,A). 

Proposition 6.1 (The area-period relation). Let A{X) be the area of the set 
{x,^ I H{x,0 < A}. Then 

(6.3) |.^(A)=T(A). 
Proof. Let w be the gradient vector field 

^^^+^^^) [-d^d-x + l^d-d'^''^ 

where p{t) = for t < | and p{t) = 1 for f > e. Then for A > e and t positive, 
exp{tw) maps the set = A onto the set H = X + t and hence 

A{X + t)= / dxd^= {exptw)* dx d^. 

JH=\+t Jh=\ 

So for t = 0, 

-A{X + t) = I dx d^ = f di{'w)dxd£, = ( L{w)dxd^. 

^ Jh<\ Jh<\ Jh=x 

But on H = X, 



d^ 
~dt 
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Hence by the Hamilton-Jacobi equations 



dH , 

ax = ^r— at 



and 



de = dt, 

ox 



the right hand side is just —dt. So 

For A = a, let c = and let 

ff^s: [0,c]->[0,a] 
be the function defined by the identities 

HUi-s) = A < 



□ 



27r 



and let 



Thus by definition 

(6.4) Acs(A) = area{iIcB<A} = A(A). 

Now let w be the Hamiltonian vector field associated with the Hamiltonian, H, as 
above and vcb the corresponding vector field for Hcb- Also as above let 7(t, A) be 
the integral curve of v on the level set, H — X, with initial point on the a;-axis, and 
let jcB{t, A) be the analogous integral curve of vcb- We will define a map of the 
set H < a onto the set Hcb < a by requiring 

i- FHcb = H, 

(6.5) ii- / maps the I'-axis into itself, 
iii. /(7(i,A)) =7cB(t,A). 

Notice that this mapping is well defined by proposition |6.1| Namely by the identity 



(6.4 1 and the area-period relation, the time it takes for the trajectory 7(i, A) to 
circumnavigate this level set H = X coincides with the time it takes for 7(73 (t. A) 
to circumnavigate the level set Hcb — X. It's also clear that the mapping defined 
by (6.5 1, i — iii, is a smooth mapping except perhaps at the origin and in fact since 
it satisfies f* Hbc = H and — vhcbj is a symplectomorphism. We claim that 
it is a C°° symplectomorphism at the origin as well. This slightly non-trivial fact 
follows from the classical Birkhoff canonical form theorem for the Taylor series of 
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/ at the origin. (The proof of which is basically just a formal power series version 
of the proof above. See |GPUj . §3 for details.) 

Now let Uq : C^(M) C°°(M) be a semi-classical Fourier integral operator 
quantizing / with the property (6.2 1. By Egorov's theorem UoS/iUq is a zeroth 
order semi-classical pseudodifferential operator with leading symbol Hqq{ ^ ) 
on the set {(a;,0 I ^qb < '^i' ^^"^ hence the operator 

is semi-classical pseudodifferential operator on this set with leading symbol of order 
h^. We'll show next that this 0{h^) can be improved to an 0{h*). To do so, however, 
we'll need the following lemma: 

Lemma 6.2. Let g be a C°° function on the set H~^{0,a). Then there exists a 
C°° function, h, on this set and a function p e C°°(0, a) such that 



(6.6) 

Proof. Let 



g = L^h + p{H). 



T(A) 



7(7(i,A)) dt 



and let gi ~ g ^ p{H) ■ Then 



T(A) 



5i(7(t,A)) dt^O. 



So one obtains a function h satisfying (6.6) by setting 



h{i{t,X))^ / gih{t,X)) dt. 



□ 



Remark. The identity (6.6) can be rewritten as 

(6.7) g^{H,h} + p{H). 

Now let —h^g be the leading symbol of 

Sn - l^lHQg{Sf^°''')Uo =: U^Rq. 

Then if h and p are the functions ( |6.6[ ) and Q is a self adjoint pseudodifferential 
operator with leading symbol h one has, by (6.7), 

^Sn-h\Ro + p{Sn)) + 0{h^). 
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Hence if we replace Uq by Ui ~ Uq expliti^Q) we have 

U.SnUl = H^QBiST) + h^P {H^^QBiSln) + 0{h'') 

= H%B{S'r) + Hl,B{S'r) + 0{h') 

microlocally on the set H^^{Q. a). 

As above there's an issue of whether (6.8 1 holds microlocally at the origin as 
well, or alternatively: whether, for the g above, the solutions h and p of (6.7) 
extend smoothly over x = =0. This, however, follows as above from known facts 
about Birkhoff canonical forms in a formal neighborhood of a critical point of the 
Hamiltonian, H. 

To summarize what we've proved above: "Quantum Birkhoff modulo fi.^" implies 
"Quantum Birkhoff modulo fi.*" . The inductive step, "Quantum Birkhoff modulo 
^2feii j]-[-^piigg "Quantum Birkhoff modulo ;i2fc+2i> jg pj-Qved in exactly the same way. 
We will omit the details. 

7. Birkhoff Canonical Forms and Spectral Measures 

Let g{s) be a function on the interval (0,cx)). Then by the Euler-Maclaurin 
formula 

Hence for / e C^(0,a), 

trace/(i^QB(5^^;^)) = / f{HQB{s,h)) ds + 0{ff-). 
Jo 



Thus by (6.11 and (6.2 1 



(7.1) i^nif) - trace/(5;,) / f{HQB{s, h)) ds + 0{h°°). 







Thus if K{t, h) is the inverse of the function Hqb{s, K) on the interval < t < a, 
i.e. for < t < a, 

K{t,h) = s ^ HQB{s,h)=t, 



then (|7.1| can be rewritten as 

= f f(f\ 

dt 



dK 

(7.2) vnU)^ / f{t)— dt + 0{hn, 



or more succinctly as 

, , dK 
(7.3) = - dt. 

Hence in view of the results of §6 this gives one an easy way to recover Hqb{s, h) 
from V and its derivatives via fractional integration. 
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8. Semiclassical Spectral Invariants for Schrodinger Operators 

WITH Magnetic Fields 

In this section we will show how the results in §3 can be extended to Schrodinger 
operators with magnetic fields. Recall that a semi-classical Schrodinger operator 
with magnetic field on M" has the form 



(8.1) 



dx 



where G C°°(M") are smooth functions defining a magnetic field B, which, in 
dimension 3 is given by i? = V x a, and in arbitrary dimension by the 2-form 
B = f^kdxk)- We will assume that the vector potential a satisfies the Coulomb 
gauge condition, 



(8.2) 



(In view of the definition of B, one can always choose such a Coulomb vector 
potential.) In this case, the Kohn-Nirenberg symbol of the operator (8.1) is given 
by 



(8.3) 

Recall that 
(8.4) 



d dp 
oxk oxk 



so the iteration formula (|2.9|) becomes 
(8.5) 



1 dhr, 



^ dt Z^ldi^dxk^''dxj'''^~' 2^ 



d dp 
' it 



dx 



dxk 



ym-2- 



from which it is easy to see that 



(8.6) 



dp dp it^ 
d£,k dxk 2 



Thus the "first" spectral invariant is 

/ Y.i^k + ak{x))^J^^\p) dxd^ = - / E ^/^P^^'^'^^ = 0, 

k k 

where we used the fact I] ff^ = 0. 



SEMICLASSICAL SPECTRAL INVARIANTS FOR SCHRODINGER OPERATORS 17 



With a little more effort we get for the next term 
4 ^ dxl 



k " 



it'' I dp dai dp dp dp d p dp ^2 

T \\^Wkd^kd^i^ ^WkWidxkdxi^ \^^d^u 

\ k^l k,l k 

—f^ dp dp dp dp 
8 ^ d^k dxk d^i dxi ' 



and, by integration by parts, the spectral invariant 

\ k k,l J 



Notice that 



and 



dxl J dxl y^dxj dxl 



So the subprincipal term is given by 



Finally Since the spectral invariants have to be gauge invariant by definition, and 
since any magnetic field has by gauge change a coulomb vector potential represen- 
tation, the integral 



is spectrally determined for an arbitrary vector potential. Thus we proved 



Theorem 8.1. For the semiclassical Schrddinger operator (8.1) with magnetic 
field B, the spectral measure ;/(/) = trace f(S^) for f £ C^(M) has an asymptotic 
expansion 

^.'"(/)^(27^;^)-"^^.;"(/);^2^ 

where 

<U)= I f{p{x,^,h))dxdi 



l//(^)(.(.,C,^))(||i.|P-2Eg). 
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9. A Inverse Result for The Schrodinger Operator with A Magnetic 

Field 



Making the change of coordinates (x,^) — > {x,^ + a{x)), the expressions (8.1) 
and (|9| simpHfy to 

'^o'if) = / ne + v)dxd^ 

and 

= f^'\e+v){\\Br-^Y.^^)dxdi. 

In other words, for all A, the integrals 

1\= I dxd£_ 

J{2+V(a;)<A 

and 

2a = / {\\Br-2Y.^)dxd^ 

are spectrally determined. 

Now assume that the dimension is 2, so that the magnetic field B is actually a 
scalar B — Bdxi A dx2- Moreover, assume that is a radially symmetric poten- 
tial well, and the magnetic field B is also radially symmetric. Introducing polar 
coordinates 

xf + X2 — s, dxi A dx2 — -ds A d9 
£,j + (,1 = t, d^i A dC2 = ^dt A 



we can rewrite the integral Ix as 

Ia = ^' / iX-V{s))ds, 



rs{X) 

2 



where V{s{\)) = A. Making the coordinate change V{s) = x s = f{x) as before, 
we get 



I;, = 7r2 t{\-x)^dx. 

dx 



A similar argument shows 



l^x = ^^ j\x-x)H{f{x))^dx, 

where 

H{s) = B{sf - AsV"{s) - 2V'{s). 
It follows that from the spectral data, we can determine 
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and 

ff(/(A))/'(A)^^^E.. 

So if we normalize V{0) = as before, we can recover V from the first equation 
and B from the second equation. 

Remark 9.1. In higher dimensions, one can show by a similar (but slightly more 
complicated) argument that V and are both spectrally determined if they are 
radially symmetric. 



Appendix A: More Spectral Invariants in 1-dimension 
For simplicity we will only consider the dimension one case. One can solve the 



equation (2.9l for the Schrodinger operator with initial conditions (2.101 and (2.11) 
inductively, and get in general 

4m 



(A.l) 

and 
(A.2) 



fe = m+l 7i + t = fc-m,,i<„ 



4 m- 2 



-'2m- 1 



fe = m+l n + t = k-m,n<,n-l 

il+-+'t=2m-l 



where a„^i and an,i are constants depending on n and h, ■ ■ • ,lt- In particular, 

/3 +4 



(A.3) 

and 
(A.4) 



(^^^^/'(x)V"(x) + ^^ev'ixr + ^^^er{x)v('Hx) + ^Q^ev^'Hx) 
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The order ti^ term is given by integrating the above formula with replaced by 
_l_ y(2;))^ By integration by parts 

j e''A{x)B^^\^ + V{x)) dxdi={-lf{2k-l)\\ J A{x)B{^ + V{x)) dxd^, 
so we can simplify the integral to 

240 + 160~ ^ 120 ^ 1440 ^ 1152 J 
Notice that 

J j-(3) ^ _ y Y'Y'" y(4) ^ y ^(4) _^ _^(5) 



and 



j V'V'V" f^^'^ = - j (2V'V'V" f''^'^ +V'V'V'V' f^^^^ 



we can finally simplify the integral to 

C2 



+ V{x)) + + V{x))j dxd^, 

J (^^{V"ix)f + ^iV'ix)^^ /(«)(| + Vix)) dxdC, 



or 

288 , 

This can also be written in a more compact form as 

e2 



It follows that 

/ ijv'v''' + %-{y'\x)f)dxd^, 

is spectrally determined. 

A similar but much more lengthy computation yields the coefficient of , which 
is given by 

^ {y'\x)f-%^{Y\x)f-%^{V\x)V\x)f-^{V\x))^f^ 



2880JV490' ' " 63' ^ " 12' ' ' ' " 144' ' " j- '2 
In other words, the integral 

is spectrally determined. 
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